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ORNSTEIN-UHLENBECK PROCESSES IN HILBERT SPACE WITH NON-GAUSSIAN 

STOCHASTIC VOLATILITY 

FRED ESPEN BENTH, BARBARA RUDIGER, AND ANDRE SUSS 


Abstract. We propose a non-Gaussian operator-valued extension of the Barndorff-Nielsen and Shephard 
stochastic volatility dynamics, defined as the square-root of an operator-valued Ornstein-Uhlenbeck process 
with Levy noise and bounded drift. We derive conditions for the positive definiteness of the Ornstein-Uhlenbeck 
process, where in particular we must restrict to operator-valued Levy processes with ’’non-decreasing paths”. 
It turns out that the volatility model allows for an explicit calculation of its characteristic function, showing 
an affine structure. We introduce another Hilbert space-valued Ornstein-Uhlenbeck process with Wiener noise 
perturbed by this class of stochastic volatility dynamics. Under a strong commutativity condition between the 
covariance operator of the Wiener process and the stochastic volatility, we can derive an analytical expres¬ 
sion for the characteristic functional of the Ornstein-Uhlenbeck process perturbed by stochastic volatility if 
the noises are independent. The case of operator-valued compound Poisson processes as driving noise in the 
volatility is discussed as a particular example of interest. We apply our results to futures prices in commodity 
markets, where we discuss our proposed stochastic volatility model in light of ambit fields. 


1. Introduction 

In this paper we introduce and analyse an Ornstein-Uhlenbeck (OU) process 

dX[t) = AX(t) dt + cr{t) dB(t) 

taking values in a separable Hilbert space H. Here, A is a densely defined unbounded operator on //, B is 
an //-valued Wiener process and a(t) is a predictable operator-valued process being integrable with respect 
to B. We shall be concerned with a particular class of stochastic volatility models rr(t) of a non-Gaussian 
nature. 

OU processes with values in Hilbert space provide a natural infinite dimensional formulation for many 
linear (parabolic) stochastic partial differential equations (see, e.g.. Da Prato and Zabczyk JH, Gawarecki 
and Mandrekar D3 and Peszat and Zabczyk li22l ~). Our main motivation for studying Hilbert space-valued 
OU processes comes from the modelling of futures prices in commodity markets, where the dynamics 
follow a class of hyperbolic stochastic partial differential equations (see Benth and Kriihner lt9l [T0l '). 

Barndorff-Nielsen and Shephard 0 proposed a flexible class of stochastic volatility (SV) models based 
on real-valued OU processes driven by a subordinator (a pure-jump Levy process with non-negative drift 
and positive jumps). This class, which we name the BNS SV model, has been applied to model financial 
time series like exchange rates and stock prices (see e.g. Barndorff-Nielsen and Shephard 0). Benth f8j 
proposed the BNS S V model in an exponential mean-reversion dynamics to model gas prices collected from 
the UK market. Later, Benth and Vos CUE) extended this to a multifactor framework to model prices 
in energy markets. Their extension of the BNS SV model to a multivariate context is based on the work 
by Barndorff-Nielsen and Stelzer 0. There are several papers dealing, both empirically and theoretically, 
with stochastic volatility in commodity prices (see e.g., Geman ED, Hikspoors and Jaimungal |[20ll and 
Schwartz and Trolle ll24l ~). 

In the present paper we lift the multivariate BNS SV model by Barndorff-Nielsen and Stelzer 0 to an 
operator-valued stochastic process, providing a very general stochastic volatility dynamics. In particular, 
we consider the ’’stochastic variance process” y(t) taking values in the space of Hilbert-Schmidt operators 
on H, 

dy(t) = cy(t) dt, + dc{t ), 


Date : June 25, 2015. 

F. E. Benth acknowledges financial support from ’’Managing Weather Risk in Energy Markets (MAWREM)”, funded by the 
Norwegian Research Council. B. Rudiger is grateful for the financial support from the Centre for Advanced Study (CAS) in Oslo. 
We are grateful for comments by Josef Teichmann. 


1 



2 


BENTH, RUDIGER, AND SUSS 


where £ is a square-integrable Levy process in the space of Hilbert-Schmidt operators on H and C a 
bounded operator on the same space. We state conditions on C and C to ensure that y is a non-negative 
definite self -adjoint operator, and in this case we define cr(t) := 3^ 1 / 2 (f). In fact, the paths of the process 
t i—^ ( C(t)f , f)n must be increasing for every / g H to have non-negative definite y. This property 
is analogous to the assumption the real-valued BNS SV model is driven by a subordinator process, since 
1 1->- (£(/)/, f)jj is equal to the scalar product of C(t) with /<g>/ in the space of Hilbert-Schmidt operators 
on H, and thus a real-valued Levy process with non-decreasing paths (i.e., a subordinator). We say that C 
has ’’non-decreasing paths” and we show that such Levy processes have a continuous martingale part with 
covariance operator having all symmetric Hilbert-Schmidt operators in its kernel. 

As a particular example a compound Poisson process is considered, where the jumps are defined to be 
the tensor product of a Hilbert space valued Gaussian random variable with itself. We demonstrate that 
such a model leads to Gamma distributed jumps for certain interesting real-valued projections of the Levy 
process. Furthermore, from a result of Fraisse and Viguier-Pla fl6l the jumps will in general be Wishart 
distributed in infinite dimensions, and we can compute the characteristic functional of C for self-adjoint 
test operators. 

Our operator-valued BNS SV model y has a convenient affine structure, and we can compute its char¬ 
acteristic function. Moreover, if C is independent of B, it is possible to derive an analytical expression for 
the characteristic function of the OU-process X(t) in terms of the semigroups associated with the drift in 
X and y and the characteristic functional of C. To achieve this result, we must impose a rather strong com¬ 
mutativity condition between the covariance operator of the Wiener noise B and the stochastic volatility 
y 1 / 2 . We find that X is affine in itself and the stochastic volatility. Also, we show that the ’’mean-reversion 
adjusted returns” of X are //-valued conditional Gaussian random variables, if these are conditioned on the 
volatility y 1 ^ 2 , which can be considered to be an observable in a simplified filtering problem (see Remark 
13.41 in Section 3). The ’’mean-reversion adjusted returns” are defined as the increments of X corrected by 
the semigroup of C. 

We relate our general analysis to commodity futures markets. In this respect, we focus on a process 
X defined on a specific Hilbert space of functions on R + , the positive real-line, and with the unbounded 
operator in the drift being A = d/dx. Then, X (t, x ) can be interpreted as the futures price at time t > 0 
for a contract delivering the commodity at time x > 0, with a dynamics specified under the Heath-Jarrow- 
Morton-Musiela (HJMM) modelling paradigm (see Heath, larrow and Morton f 19| and Musiela (22). We 
connect our general S V modelling approach to the analysis in Benth and Kriihner 19ll 1 011 and the ambit field 
approach in Barndorff-Nielsen, Benth and Veraart (3l[4l. We remark that this discussion can be extended to 
forward rate modelling under the HJM paradigm in fixed-income theory (see Filipovic El and Carmona 
and Theranchi fl3l for an analysis of HJM models in infinite dimensions for fixed-income markets.). 

Our results are presented as follows: In the next section we introduce the operator-valued BNS SV 
model and analyse its properties. Section 3 defines the volatility-modulated OU process X along with a 
discussion of its characteristics. Finally, in Section 4, we discuss our model X in the context of commodity 
futures price modelling. 

2. Operator-valued BNS stochastic volatility model 

Throughout the paper, (0, T. {Pt\t> o, P) is a given filtered probability space. Let H be a separable 
Hilbert space with inner product denoted by ('I ')h and associated norm | • //. Introduce 'H := L H s(H), 
the space of Hilbert-Schmidt operators on H into itself, with the usual inner product denoted by (■> ’)« and 
associated norm || - \\u- As H is a separable Hilbert space, 'H becomes a separable Hilbert space as well. 

Introduce C G IX'H), that is, a bounded linear operator from 'H into itself. In this paper, we shall pay 
particular attention to two specific cases of C, namely, the operator 

(2.1) Ci T^CTC* 

or the operator 

(2.2) c 2 -.n^n, r^cr+rc*. 

Here, C € L(iT), L(H) denoting the space of bounded linear operators in H into itself. We shall exclu¬ 
sively focus on C ^ 0. The following lemma provides us with crucial properties for C,;. i = 1,2: 
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Lemma 2.1. It holds that (Cj G L {PL) for Cj defined in ( 12.1b and (12.2b . with |[Ci|U, < \\C\\l p and\\C. 2 \\op < 
2||C|| op . Moreover, (C fT)* = CfT* for every T GPL and i = 1,2. 

Proof. For S,T GPL, C fiS + T) = C,*S + C,T, where i = 1,2. Hence, linearity holds. Moreover, for an 
orthonormal basis {e„} ra gN in H, 

HCrTI^ = \\CTC*\\ 2 n 

oo 

= J2 \CTC*e n \ 2 H 

n =1 

oo 

n =1 
oo 

<||C||4 p ^|re„| 2 ff 

n=l 

= iicii 0 4 P imi^. 

Here we have used that ||TC*||^ = ||CT*||«. For C 2 , we have by the triangle inequality, 

\\C 2 T\\n = \\CT+TC*\\ n 

< \\CT\\n + \\TC*\\n 
<l|C||op||r|k + ||C|| op ||Tl« 

= 2||C||op||7l«. 

Hence, the first claim of the lemma holds. 

For T G PL , it follows that 

(C 1 Tf,g) H = (CTC*f,g) H 
= {f,CT*C*g) H 
= (f,C 1 T*g) H . 

An analogous computation shows that also CfT = C 2 T*, and the second claim of the lemma holds. 
Hence, the proof is complete. □ 


Since C G L{PL), it follows that C generates a uniformly continuous Co-semigroup §(f), t > 0, with 
Ei(t) = exp(fC) (see for example Gawarecki and Mandrekar lfl7l Thm. 1.1]). We note the following for 

C i,i = 1,2: 

Lemma 2.2. For the Co-semigroup S i generated by C, in (12.1b and (12.2b . i = 1,2, resp., we have 

00 4 .71 

§ 1 (^)r=y — c n T(c*r, 

n\ 

n —0 

and 

S 2 (t)T = exp {tC)T exp(fC*), 

for every T G PL- 

Proof For T GPL, we find for n > 1 

CIT = Cl~ l (CTC*), 

and iterating this yields 

C"T = C n TC* n . 

Hence, the result for §1 follows. 

For the case C 2 , note that 

00 4 .n-\-m 

exp(fC)T exp(fC*) = V —-C n TC* m . 

^ n\m\ 

n,m —0 
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On the other hand, 


oo j-k 

ex P (iC 2 )T = ^ n QT. 

k =0 

Spelling out T and comparing with the terms in the double-sum above, we show the second result, 
proof of the lemma is complete. 


The 

□ 


We now introduce the operator-valued BNS stochastic volatility model. To this end, assume that 
{4’(f)}t>o is a -valued stochastic process satisfying the dynamics 

(2.3) dy(t) = c y(t) dt + dc(t) t(o) = y 0 . 

Here, C is an %■-valued Levy process and To G PL. We suppose that C is square-integrable, with covariance 
operator Q^. Recall that Qc is a self-adjoint non-negative definite trace class operator on PL. We have. 

Lemma 2.3. For every t > 0, it holds 

f* m - snTwL («) ds < ^^(e 2t||c|1 - -!) < 00 • 

Proof. Note first that for any T £ PL. we have by the representation of S, 

l|S(«)T|| w < ||S(u)||op||T||w 
00 k 

simi^Elr wt 

k =0 

= e“H c ll° p ||T||-H • 

But then, for an orthonormal basis {T n }neN C 'H , 

OO 

= E 

n= 1 

00 

<e 2 “ll c H<*^||Qj/ 2 7 ;||^ 

n—1 

= e 2 “ l|c||o »Tr(Q z; ). 

Here we have used the fact that 

OO OO 

Tr(Q £ ) = ^(Q c T n ,T„)n = E \\®c 2 T*\\h ■ 

n= 1 n=1 

Hence, since ||C|| op < 00 and is a trace class operator, the result follows. □ 


Invoking this lemma, it follows from the theory of Hilbert-space valued stochastic differential equations 
(see e.g. Peszat and Zabczyk | [22l ) that there exists a unique mild solution to (12.3b 

(2.4) y{t) = §(t)y 0 + f §{t - s) d£(a), 

Jo 

for t > 0. In the next lemma we derive a bound for the L 2 - norm of y: 

Lemma 2.4. It holds that 

E[||T(f)||«] < ce 2tn ” 

for a constant c which is bounded by max(2||To ||^, Tr(Q £ )/||C|| 0/; ). 

Proof. From the mild solution of y(t) in (12.4b and the triangle inequality we find, 

E[||T(i)|&] <2||S(t)To|& + 2 l* m-s)®]/ 2 \\ 2 LHs{H) d s , 

Jo 

where we used Cor. 8.17 in Peszat and Zabczyk l22l . But from Lemma [2731 above the result follows. □ 
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Let us compute the conditional characteristic function of y(t): To this end, let t > s and note that y(t) 
given ;y(s) has the representation 

(2.5) y(t) = S(t — s)y(s) + f 8(7 — u) d,C(u). 

J s 

Before proceeding, we recall the cumulant of £, that is, the characteristic exponent of the Levy process £ 
defined to be E[exp(i(£(7), T)n\ = exp(7\D,c(7”)) for T £ 77 (see Peszat and Zabczyk ll22l Thm. 4.27]): 

(2.6) *c(T) = i(V 1 T) n -\{<QPcT,T) H + [ -\-\\\\z\w<x{Z,T)n)y{dZ). 

J 'LL 


Here, following Peszat and Zabczyk [;22] Thms. 4.44 and 4.47], v is the Levy measure on 77 and V £ His 
the drift of the Levy process, where for T € 77, 

E[(£(l), T)u] = (V, T) u + [ (2, T)u v(dZ). 

Ju\n\z\\-u<n 


Furthermore, and 


(QcT,U)h= [ (T,Z) n (U,Z) n v(dZ), TM&H. 

Jn 

We have the following proposition, showing that y is an affine process in 77,: 

Proposition 2.5. For any T € H it holds that 

E UTMfOw | j-1 = exp fi(3^(s),S*(f — s)T)u + f ^ c ^*{u)T)du) . 


Proof. From (12.5b we find for T € 77, 


E 




= e i<s(*-s)T(s),r>KE 
= e i<T(^),S*(i-s)T)„E 


s(t-u) dC(u),T)-H I jr 


d.C(u),T)n 


Here, we have appealed to the independent increment property of Levy processes. Hence, from Peszat and 
Zabczyk l22l Thm. 4.27] it holds that 

pt — S 


E 




= exp 


'Tc(S*(u)T) ds 


with ^c defined in (12.6b . The result follows. 


□ 


To define a stochastic volatility based on y in ( 12.4b we must impose positivity constraints. This means 
that we want to restrict our attention to 3^’s which are self-adjoint, non-negative definite Hilbert-Schmidt 
operators on H for each t > 0. We now analyse additional conditions on C and £ ensuring non-negative 
definiteness of y. First, we show that y[t) is self-adjoint whenever £(f) is under a mild condition on C: 

Proposition 2.6. Suppose that (CT)* = CT* for any T £ 77. If {£(7)}oo is a family of self-adjoint 
operators on H and 3^o L self-adjoint, then y(t) is a self-adjoint operator on H for every t > 0. 

Proof Let f,g £ H. Then we compute, using the dynamics of y in (12.3b . the assumption on C, the 
self-adjointness of £(f) and the definition of Bochner integration: 

( y(t)f,g) H = f (c y(s)f,g) H ds + {C(t)f,g) H 

Jo 

(f,Cy*(s)g) H ds + (f,£(t)g)H. 

Thus, as /, g £ H are arbitrary, we find that 

dy*(t) = Cy*(t)dt + d£(t), 
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with initial condition y* (0) = Vo • But by uniqueness of solutions of this linear stochastic differential 
equation, y*(t) =y(t). □ 


Recall from Lemma lZTl that (C,;T)* = C fT* for * = 1,2. 

Example 2.7. A trivial way to introduce a self-adjoint Levy process C in FL is to take any real-valued Levy 
process L and multiply it with a self-adjoint operator IA £ 7 j, i.e., C(t) = L(t)U. For S, T £ FI, 

E [<£(i),S)„<£(i),T>«] = E [L 2 (t)}{U,S) H {U,T)n = tVar(L(l))(U® 2 S,Tj n . 

Thus, the covariance operator for this Levy process becomes Qc = LA® 2 , i.e., the tensor product of U with 
itself. We show that (Q>£ is a self-adjoint, non-negative definite trace class operator. Indeed, it is obviously 
linear and 

\\QcT\\n = mil,r)n\\n<\mn\\T\\n, 

which shows Q c £ L(FL). Moreover, 

(QcS,T) n = (U,S)h{U,T)h = (S,LA® 2 T)u = (S,Q c T)n 

and 

(QcS,S) n = (U,S) 2 l >0, 

which show that Qc is a self-adjoint and non-negative definite operator on FL. Finally, for an orthonormal 
basis {7n}neN in FL, 

oo oo 

Tr(Q £ ) = J2^cT n ,T n )n = X>,T n & = Mn 

n =1 n =1 

where we used Parseval’s identity. Hence, Q c is trace class. Of course, if we add the assumption thatU is 
positive definite and L(fi) is taking values on K+,0 it follows that 

(C(t)fJ) H = L(t)(UfJ) H >0, 

for any f £ H, and thus C{t) is non-negative definite. 

This simple example of an operator-valued Levy process C brings us to the question of non-negative 
definiteness of y, which we investigate next. First, let us define what we mean by non-decreasing paths of 

C: 

Definition 2.8. We say that the FL-valued Levy process C has non-decreasing paths if C{t) is self-adjoint 
and t i—>■ ( C[t)f , f)n is non-decreasing in t > 0 for every f £ H, a.s. 

Note that as £(0) = 0 by definition of the Levy process, the non-decreasing paths property implies 
> 0 for every t > 0, a.s.. But then it follows that C(t) is a non-negative definite operator. In 
fact, something slightly stronger holds: 

Lemma 2.9. Assume C is an FL-valued Levy process with non-decreasing paths. Then C(t) — C(s) is a.s 
non-negative definite for every t > s > 0. 

Proof. For t > s > 0, we have for / £ H 

am - C(s))f , f) H = (£(*)/, f) H - (£(*)/, f)H 

which is non-negative a.s by the non-decreasing path property of t i-» ( C(t)f , f)n- The assertion fol¬ 
lows. □ 


As we shall see, this monotonicity property of the ’’paths” of C is exactly what we need in order to show 
that y(t) is non-negative definite for every t > 0. But first, let us do some analysis of Levy processes in FL 
with non-decreasing paths. 


l 


This means that the Levy process is a so-called subordinator, that is, a process with only positive jumps and non-negative drift. 
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Define for the moment Lf(t) := ( C(t)f , /7 for given f £ H. We show that this is a Levy process 
on the real line. To this end, consider the functional J 7 / : 7i —> R defined as T/(T) = (T/, f)n- This is 
obviously linear, and since 

|T/(T)| = 1(77,/) ff | < IT/71/7 < IITU/7, 

we have T/ C PL* . Hence, there exists a unique element in 7f, which we also denote by T/, 

T/(T) = {T,T f ) n . 

In the following we identify Hilbert-Schmidt operators on H with H (g) H. Similarly, the Hilbert-Schmidt 
operator h* h for h £ H is written as h 0 h. Then for any Hilbert-Schmidt operator V we have the 
following identity 

(V,h®h) H = (Vh,h) H - 

Thus, we have J 7 / = / ® /. Indeed, a straightforward calculation shows, 

OO 

11/ ® f\\h = ® f)e n , (/ ® f)e n )H 

n= 1 
oo 

= X]((/’ e ")ff/’ (/’ e n ) Hf)H 
n= 1 

oo 

= I/Iff y~l(/> e n)g 

n=l 

= l/ltr ■ 

Hence, f <8> f £ H with norm ||/ <g> /|7 = I/I 2 h- Furthermore, 

OO 

(T, f®f)u = ^^(Te„ , (/ <8> f)e n ) H 

n— 1 
oo 

= y~l(Te„, (/, e n ) H f)H 

n— 1 
oo 

= y ~1(T(/, e n )ge n , f)n 

n= 1 

= (T/,/7- 

By definition of an 7L-valued Levy process, f i—>• (C(t), T)n is a real-valued Levy process for any T £'H. 
Therefore, in particular, Lf(t) = (. C(t)f, f)n is a real-valued Levy process by choosing T = / <S> / . If, 
furthermore, £ has non-decreasing paths, it follows that L/ is a Levy process with non-decreasing paths, 
i.e., a subordinator. We have the property of the continuous martingale part of C: 

Proposition 2.10. Let C be the Levy process defined after (12.3b with non-decreasing paths, and denote 
the covariance operator of the continuous martingale part by Q^. Let T be a symmetric Hilbert-Schmidt 
operator. Then Q^T = 0, that is, T £ ker(Q^). 

Proof. Let first T = f ® f with f £ H. Then the continuous martingale part of the characteristic function 
of L f (t) = (C(t),T f ) n is <Q£T>,T» w , which must be zero due to the non-decreasing paths of Lf(t). 
But then 

(Q!cT>,.t»« = ||(Q!c) 1/2 t>||^ = o, 

and thus T/ is in the kernel of (Q 1 ^) 1 / 2 . As it holds, 

Qc^f = (Qc) 1/2 (Qc) 1/2 ^f = (Q£) 1/2 o = 0, 

we can conclude that J 7 / £ ker(Q>£). 

Now let T be a symmetric Hilbert-Schmidt operator as in the proposition. It can be shown that T must 
be of the form 

T = ^2 7 k,iek ® ej, 

k,leN 
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with Y2k 1 7k i < 00 an d 7k,i = 7 i,k, see Lemma lATI for a sketch of the arguments. Therefore we can write 

"T — ^ ^ 7fc,fc^fc ® &k H - ^ EE 7fe,;(efc <g» ez + e; <S> efc) 

*eN AeN l<k 

= E^®e t +2EE 7 m (( e fc + e/) ® (efc (8> e{) - e k ® e k - eo 8> ez). 

fcGN fcGN i<fc 

With this we compute 

Q/:T = ^2 7k,kQc( e k ® e fc ) 

fcGN 

+^EE 7k,i(Q°c(( e k + ez) (8) (efc + ej)) - Qj(efc <8> e fc ) - Q^(ez (8) ez)), 

fcGN Z<fc 

which ends the proof since applied to a symmetric operator is zero by the first part of the proof. □ 

As the space of symmetric Hilbert-Schmidt operators does not span PL, we cannot conclude that = 0, 
i.e., that C does not have a continuous martingale part. Recall that subordinators on R do not have any 
continuous martingale part. 

Denote now by PL+ the convex cone of non-negative definite operators in 'H. 

Proposition 2.11. Assume that C(PL+) C PL+. If C(t) is an PL-valued Levy process with non-decreasing 
paths and is non-negative definite, then y is non-negative definite. 

Proof Recall that 

y(t) = s(t)y 0 + [ s(£ — s) c(s). 

Jo 

It holds, 

OO 

sm = e tc y Q = y jfi£ k yo, 

fc=o 

which is then a non-negative definite operator by the assumption on C. 

Next, we know that f* S(t — s ) dC(s) is defined as the strong limit of XE=i — s m) A£(s m ) in 
L 2 (fl x [0, £];?£). Here, {s m }m=i is a nested partition of [0, f], and A C(s m ) ■= £(sm+ 1 ) — £(s m ) is 
an increment of C. But A C(s m ) is non-negative definite a.s. by Lemma [2~9l and therefore each term in 
the sum above is positive, a.s., since § preserves non-negative definiteness by assumption on C. Hence it 
follows that f ( j S(t — s) dC(s) is non-negative definite a.s., and the proof is complete. □ 


From Pronosition l2.6l and Proposition 12.1 H it follows that under the assumptions 

a) (cry = c r*, 

b) C (H+) C U+, and 

c) C(t) is a self-adjoint and non-negative definite square-integrable Levy process with values in PL, 

then y(t) becomes a self-adjoint, non-negative definite square integrable process with values in PL as long 
as the initial condition A'o is self-adjoint and non-negative definite. Hence, we have a unique square root, 
y x l 2 ( t ) for every t > 0. We shall use this to model the stochastic volatility. 

Lemma 2.12. It holds that C±(Pl+) C PL+. 

Proof. We recall the definition of Ci in (12.l) . Let T G PL+. Then, for any / G H 

(c iTf,f) H = ( crc*fj) H = ( rc*f,c*f) H > o. 


Hence, the result follows. 


□ 
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In fact, for C 2 we cannot prove that it preserves the property of non-negative definiteness. But recalling 
the proof of Prop. l2.in it is indeed the associated semigroup of C that must be non-negative definite. As we 
have that §2 (i)T = exp(/C)Texp(/C*) from Lemma [Z2l it follows that § 2 (/)(?/+) C "//+, and we can 
conclude that y with C = C 2 is also non-negative definite whenever C has non-decreasing paths and To is 
non-negative definite. Indeed, by inspection of the proof of Prop. 12. in we can substitute the condition b) 
<C(7/+) C 7/+ on C with the condition 
h’) CH+, t> 0. 

In conclusion, if we use C = Cj for either i = 1 or i = 2 in the definition of the volatility process y, 
we obtain a non-negative definite operator under appropriate conditions on C and To- We recall that the 
choice C = C 2 can be seen as the analogue of the matrix-valued volatility model by Barndorff-Nielsen and 
Stelzer Q. 

Let us discuss the particular case when C is a compound Poisson process. To this end, we define 

N{t) 

m = Y, X >-' 


(2.7) 


*=1 


where TV is a real-valued Poisson process with intensity A > 0 and are i.i.d. square-integrable 

//-valued random variables. Note that for f £ H, we find from the linearity of the inner product 

N (t) N (t) N (t) 

(m,f 0 f) n = <E f ® /)« = E / ® f)n = E (^/. f)n • 

i=l Z =1 2=1 

Hence, L/(f) := f <g> /)^ is a real-valued compound Poisson process with jumps given by the i.i.d 

random variables {Xif, f ) n ■ The process Lf(t) has non-decreasing paths if and only if X is self-adjoint 
and (X /, /) n is distributed on R + , where the latter holds if and only if X is non-negative definite, i.e., 
X £ 'H + . Next, introduce the map 4>f : T~L + —> R + by 

MZ) = (ZJ®f)n- 

For any Borel set A C R+, we define P t j )f (A) := Px{4 >^ 1 (A)) where Px is the law of X. But then 


P uz P^ f {dz) = 


[ (e iu ' o cf> f )(Z)Px{dZ) = [ e iu < z ’ f ®»*Px(dZ), 

JH+ JH+ 


and P Xf ( A ) = P x (<p f 1 (A)) with P Xf being the law of X f := (Xf, f) H . 

Suppose that Z is an //-valued centred square-integrable Gaussian random variable with covariance 
operator Qz- Let Xi = Zf 1 . i = 1,2,...,, where { Z ,} are independent copies of Z. First, it is simple 
to see Z® 2 is also a Hilbert-Schmidt operator, that is Z® 2 £ H, since 

OO OO 

l|0 S2 ||w = ^ \Z® 2 e n \ 2 H = c^hIZIh = \Z\% < 00 . 

n =1 n =1 

This '//-valued random variable has expected ('//-valued) value E[Z® 2 ] = Qz, which can be seen from the 
following calculation: given T £ 7~L, then by linearity of the expectation operator 

OO 

<E [^® 2 ] ,T)u = ^(E[Z® 2 ]e n ,re„) ff 

n= 1 
00 

= [(Z® 2 e n ,Te n ) H \ 

n= 1 
00 

= ^2 E [(Z, e n )H(Z, Te n )rr] 


n—1 

00 


= y^(Qze n , Te n )n 


n—1 


= (Qz,T)h. 
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Furthermore, Z® 2 is self-adjoint and non-negative definite, since (Z® 2 /, /)# = (Z, f) 2 > 0. From this 
we also see that the jumps of Lf(t), the compound Poisson process £ evaluated at / <g> /, is given by 
(ZJ)jj, with ( Z,f)u being a real valued centred Gaussian variable with variance \Q^ 2 f\%. Hence, 
(ZJ) 2 h becomes Gamma distributed with scale parameter 2\Q 1 J 2 f\ 2 H and shape parameter 1/2. In fact, 
something much more general can be said about the compound Poisson process £ for jumps given by 
X = Z® 2 . Indeed, if T € 7/ is self-adjoint, then it follows from Prop. 3 in Fraisse and Viguier-Pla m 
that the characteristic functional of ( Z® 2 , T)u is, 

(2.8) E [exp(i(Z 02 , T)u)\ = (det(J - 2i TQ Z ))~ 1/2 ■ 

Here, I is the identity operator on H and det is the Fredholm determinant. We can interpret Z® 2 as being 
infinite dimensional Wishart distributed. By conditioning of N(t) and appealing to the independence of 
the jumps Xi, we find the cumulant \p£ of £ defined in (12.6b to be 

(2.9) * C (T) = InE [exp(i(£(l), T)h)\ = A ((det(J - 2iTQz))~ 1/2 - l) , 
for any self-adjoint T £ PL. 

Suppose now more in general that Z is an //-valued centred square-integrable random variable. Then Z 
has a self-adjoint non-negative definite continuous linear covariance operator Qz, too. Let X, = Z® 2 , i = 
1,2,...,, where {Zi }i £ N are independent copies of Z. Then by the same calculations as before ||Z <82 ||-h = 
\Z\jj < oo and (E [Z® 2 ] ,T)u = ( Qz,'T)n ■ Also here Z® 2 is self-adjoint and non-negative definite, 
since = (Z, f) 2 > 0 and the jumps of £/(/), the compound Poisson process £ evaluated at 

/ < 8 > /, is given by ( Z,f ) 2 H , with (Z, f)n being a real valued centred variable with variance 1 0%“ f \ 2 H . The 
cumulant has then to be computed for each case separately. 


3. A VOLATILITY-MODULATED ORNSTEIN-UHLENBECK PROCESS 


Let X be a stochastic process with values in H given by the Ornstein-Uhlenbeck process 
(3.1) dX(t) = AX(t)dt + y 1 / 2 (t)dB(t) X(0)=X o . 

Here, B is a //-valued Wiener process with covariance operator Q, which is a self-adjoint, non-negative 
definite trace class operator on H. Furthermore, Xq £ H and y is given in (12.4b . being the solution of the 
dynamics ( 12.3b from the previous section, where we assume that To is self-adjoint, non-negative definite 
and £ is a 7/-valued Levy process with non-decreasing paths. We suppose that (C T)* = C T* for every 
T £ W and C (TL+) C %+ (or, that the semigroup §(/) of C has this property). Then by Props. [2~6l and 
12.111 y(t) is self-adjoint, non-negative definite, and we can define its square root T 1 ^ 2 (/)- Finally, A is a 
(possibly unbounded) linear operator on H, densely defined, generating a Co-semigroup S. 

Let us first show that the stochastic integral in (13.1b makes sense. The following proposition is crucial: 


Proposition 3.1. For every t > 0, it holds that 


E 


Tr(Q 1 / 2 T(/)Q 1/2 )l =Tr(Q 1 / 2 S(/)ToQ 1/2 ) + Tr(Q 1 / 2 / S(a)daE[£(l)]Q 1 / 2 ) 


where f* S(s ) ds is the Bochner integral of s §(s) £ Lhs(H) an d E[£(l)] is the operator-valued 
expected value o/£(l). 

Proof. First, note that the trace is linear, to give 

E [Tr(Q 1 / 2 T(/)Q 1/2 )l =Tr(Q 1 / 2 §(/)ToQ 1/2 )+E Ti/Q 1 / 2 f S(i - s) dC{s)Q 1/2 ) 

L J L Jo 

Suppose for a moment that X is a '//-valued integrable random variable. Then 


E 


Tr(Q 1 / 2 A'Q 1 / 2 )] =J2 K [(Q 1/2 XQ 1/2 e n ,en) H 

n —1 
oo 

= Y J A( X Q 1,2e n,Q 1/2 e n )H 
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But (Xf,f) H = {X , / ® f)ui which holds due to the isometry of the Hilbert-Schmidt operators with 
tensor products of Hilbert spaces, and 


mxf, f) H ) = E[{X, / ® />*] ={M,f <8 f)n , 

for some M £ H. This operator is called the mean of X, and we write E[X] = M, the operator-valued 
expectation of X. Thus 


E 


Ti{Q 1/2 XQ 1/2 ) 


= 'STe[(X, Q 1/2 e n ®Q 1/2 e n )-H 

n= 1 
oo 

= Q l/ 2 e n ) H 


n—1 


= ^(AfQ 1/2 e„,Q 1 / 2 e n ) 

n= 1 

= Tr(Q 1 ^ 2 MQ 1/2 ) 

= Tr(Q 1/2 E[-T]Q 1/2 ). 


Letting X = f* S(t — s ) dC(s), we hence obtain 

E[Tr (Q 1/2 - s) dC{s)Q l/2 )\ = Tr ^Q 1/2 E J S(t-s)dC(a) 

We derive an expression for the mean of the stochastic integral. 

Recalling (the proof of) Prop. 12.51 we find that with 6 £ R 


Q}! 2 


E 


e i(/o S(t-s) d.C(s),6T)u 


= exp ( / ^£(S*(u)(0T)) du 


= exp (/ 'p£(0§*(u)(T)) du 

with defined in (12.6b . Since ’f',c(0) = 0, we find 


d ^ 
de E 


j(fo S(t—s) d£(s),0T>n 


|e=o= / ^^(«S*(«)T)du|9=o. 


But, for any S £ H, 

^c{0S) = i(V, S)n - e(Q° c S,S) H + i J «Z- <Z, S)«1(||Z|| W < 1)) t/(dZ). 

Therefore, 

(f S(i - s) dC{s ), T)« = (—i) J* ^$ £ (flS*(«)T)k=o du 

= f (P,S*(u)T)«+ f (Z, §* {u)T)n v{dZ) du 

Jo J \\Z\\w>l 


E 


(u)27,r>« 


(u)Z, T)n v{dZ) du 


ll-Z||w>l 


= (/s (u)(V+ f Zv(dZ)du),T)u 
Jo J||^||tk >1 

= {[ §{u)du(V+ [ Z v(dZ), T)n ■ 
Jo J\\Z\\ V >1 


Thus, since E[£(l)] = 12 + J|| 2 || M>1 Z v{dZ), we get 

t 


E 


§(f — s) d£(s) 


Wo 


= [ S(u) duE[£(l)]. 
Jo 
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This completes the proof. 


□ 


To have the stochastic integral f* y 1 ^ 2 (s) dB(s) well-defined, the integrand must satisfy the condition 


(3.2) 


E 



\\y 1/2 (s)Q 1/2 fnds 


< oo. 


But \\y 1 / 2 {s)Q 1 l 2 \\1 l = Tr(Q 1 / 2 J^(s)Q 1 / 2 ). From Prop . 13.11 above we see that the expected value of this 
trace is integrable in time on any compact set. Thus, y 1 ' 2 can be used as a stochastic volatility operator in 
the dynamics of X in (13.1b . 

If the stochastic integral f* S(t — s)3^ 1 ( 2 (s) dB(s) exists, then we have the mild solution of ( 13.1b 

(3.3) X{t) = S(t)X 0 + f S(t- s)y 1/2 (s)dB{s), 

J 0 

for a given initial condition X(0) = Xq E H. The stochastic integral is well-defined since 


m - s)y^ 2 (s)Q 1 / 2 \\ n < ||S(i - s)||op||3 ;1/2 (s)Q 1/2 ||w . 

By Yosida l26l , the operator norm of the semigroup S is at most exponentially growing. Hence, in view of 
Prop. 13. ll integrability holds. 

Here is a result on the characteristic function of the process X(t): 


Proposition 3.2. Suppose that there exists a self-adjoint, positive definite operator T> £ /.(//) such that 
y 1/2 {s)Qy 1/2 {s) = V 1 / 2 y(s)V 1 / 2 for all s > 0. Then, if C is independent of B, 


E 




= exp (i(X 0 ,S - 2 CVo, J S *(s)((V^ 2 S*(t - s)f) <g> (V 1 / 2 S* (t - a)/)) ds) H 

x exp (&£ i J ^*( s — u )(V 1 ^ 2 S*(u)f ® V 1 ^ 2 S*{u)f) dii'j ds'j , 


for any f £ H. 


Proof First, from the mild solution of X (t) we find for / £ // 


(X(t), f) H = ( S(t)X 0 ,f) H + ( r s(t - s)3 ;1/2 (s) dB(s), f) H . 

Jo 

We compute the characteristic function of the random variable (f* S(t — s)y 1 ^ 2 (s) dB(s ), /)h- Since £ 
and B are independent, we have that y and B are independent. From the tower property of conditional 
expectation, we therefore get after conditioning on the er-algebra generated by the paths of y-. 


E 


6XP (^/ S<(L ~ s X /2 ( s ) dB ( s )J) H 


= E 


exp {Qy 1/2 (s)S*(t- s)f,y 1/2 (s)S*(t- s)f) H ds 


From the property of the operator V, 

( Qy 1/2 s*(t - s )f,y 1/2 (s)s*{t - 8 )f) H = {y 1/2 {s)Qy 1 ' 2 {s)s*(t - s)f,s*(t - 8 )f) H 

= {v 1,2 y{ s )v 1/2 s*{t - s)f, s*(t - 8 )f) H 
= (y{s)v 1/2 s*(t - s)f,v 1 / 2 s*(t - 8 )f) H 
= (y(s), (© 1/2 5*(t - s)f ) ® {VV 2 S*{t - S)f))u . 


For simplicity, introduce for the moment the notation T(s) £ ‘H for the family of operators parametrized 
by time s > 0, defined by 

T(a) = (V^S^s)/) ® (DV 2 5*(a)/). 

Thus, from the mild solution of y. 
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We have that 


f (S(s)y 0 ,T(t-s)) n ds = (y 0 , f §*( S )T(t~s)ds)u 
Jo Jo 


where the integral on the right-hand side is interpreted in the Bochner sense. It holds, after appealing to a 
Fubini theorem for stochastic integrals in Hilbert space (see Peszat and Zabczyk (I 22i Theorem 8.14]) 

f (f §(s — it) dC(u),T(t — s))n ds = f (f S*(s — u)T(t — s) ds, dC{u))u . 

Jo Jo Jo Ju 

The ds-integral inside the inner product is again interpreted as a Bochner integral. Hence, 

r / i r* 

E 


exp J (J S {s-u)dC{u),T(t-s))uds 

GXP / (/ ~ u )T{t~ s)ds,dC( u ))u 


= E 


= exp 


= exp 


'i'c - 


4 'c - 


-^J S*{s~u)(V 1 / 2 S*{t-s).f^V 1/2 S*(t-s)f)ds^j du) 
§*(s-u)(T» 1/2 5*(u)/®X> 1/2 5*(it)/)dM^ ds ^ . 


This proves the Proposition. 


□ 


The result shows that we recover an affine structure of X in terms of Xq and To ■ Note that if Q 
commutes with T(s), then Q 1 / 2 commutes with T 1 ^ 2 (s), and we find 

y 1/2 {s)Qy 1/2 ( s ) = q 1/2 t(s)q 1/2 . 

Hence, in this case T> = Q. Indeed, this puts rather strong restrictions on the volatility model y. A 
sufficient condition for y commuting with Q is that Q commutes with To an d C{t) for all t > 0, and that 
C(T) Q = C(7” Q ) and QC(T) = C(QT) for every T £ Ji. If this is the case, we have from the dynamics 
of V in (O 

(3.4) Qy(t) = Qy 0 + [ C(ST(s)) ds + QC{t ), 


and 

(3.5) 

Introduce now the notation 


y(t)Q = QTo + [ C(y(s)Q)ds + QC(t). 
Jo 


(3.6) 


C Q {t) := QC(t), 


which is an 'H-valued process. It is in fact a Levy process with values in H. Indeed, its conditional 
characteristic function is (here T £ T~L and t > s ) 


E 


3 i(£ a (t)-£ c (s),T>w | 


= E 
= E 
= E 


e HQ(C(t)-C(s)),T)-H | jr. 
e i (C(t)-C(s),QT)-H | J7 
J (C(t)-C(s),QT)n 


= exp ((t - s)'Sc(QT)) 


by the independent increment property and the definition of the cumulant of C. Hence, C.Q.{t) — Cq(s) is 
independent of F s with a stationary distribution, which implies that Cq is a Levy process. Its covariance 
operator is given by Q,c q = QQrS, which is easily seen from 

E [{C Q (t),T)n(C Q (t),S)n} = E[<£(i), QT)«(£(i), QS) n ] 

= (Q cQT, QS) U 
= (QQcQT, S) n , 
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with T,S £ R. Therefore, we have a mild solution of the equation for jpg := QjV in ( 13.4b given as 

(3.7) y Q (t) = S(t)Qy 0 + [ S(t - s ) dC Q (s). 

Jo 

Moreover, we see that y(t) Q in (13.5b solves the same equation, and thus Qy(t) = y(t) Q by uniqueness of 
solutions, and the claimed commutativity follows. We remark that if Q commutes with C, then the assumed 
property of C = C, holds for i = 1,2. Also, if C is the simple choice as in Ex. 12.71 it commutes with Q 
whenever U commutes with Q. 

Let us investigate the ’’adjusted returns” implied by the model. To this end, fix At > 0, and define the 
’’adjusted return” by 

R(t, At) = X(t + At) - S(At)X(t). 

From ( 13.3b , we find after using the semigroup property of S, 

/ t+At 

S(t + At-s)y 1/2 (s) dB(s). 

We have: 

Lemma 3.3. Let be the a-algebra generated by the paths ofy. Then R(t , AtjlJ 7 ^ is a mean zero 
H-valued Gaussian random variable, with covariance operator 

S(t + At- s)y 1 / 2 (s)Qy 1 / 2 (s)S*(t + A t-s)ds. 

Proof. By inspection of the proof of Prop. 13.21 we find for f € H 

( i /'£+A£ 

~2 J \Q 1/2 y 1/2 S*(t +At - s)f\ 2 H ds 

This is the characteristic function of a Gaussian mean-zero real valued random variable. Hence, f?(f, At)\J- y 
is Gaussian in H with mean equal to zero. The conditional covariance operator follows by a direct compu¬ 
tation. □ 



The stochastic volatility model yields a Gaussian variance-mixture model for the adjusted returns (see 
Barndorff-Nielsen and Shephard J6] for mean-variance mixture models and stochastic volatility in finance). 
Remark that if Q and y commute, the conditional covariance operator becomes 

S(t + At — s)y Q ( s)S* (t + At — s) ds , 

with the definition of jpg given above. 

Remark 3.4. In Lemma 13. jj a simplified filtering problem for the adjusted returns of the model is solved. 
Here the observable is the volatility y 1 / 2 . Compared to more general filtering models, as e.g. those 
described in Xiong 11251 . our filtering model is simple, as the observable does not depend on the signal X, 
for which the adjusted returns are computed. 

4. Application to forward price modelling 

Let H = H w , the Filipovic space of all absolutely continuous functions / : R + -+ R such that 

/»00 

(4.1) \f\l ■-f(0) 2 + w(x)\f'(x)\ 2 dx < oo , 

Jo 

where w : R+ —>• M+ is an increasing function with w( 0) = 1. We assume that J 0 °° w~ 1 (x) dx < oo, 
and denote the (naturally defined) inner product (•, •),,,. It turns out that H w is a separable Hilbert space 
equipped with the norm | • Moreover, the evaluation functional 6 x (f) = f(x) is continuous on H w . As 
a linear functional, we can express 8 X by (•, h x ) w , with 

rxAy 


(4.2) 


h x (y) = 1 + 


0 


w 1 ( 2 ) dz, y £ R + . 
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See Filipovic |fl5l for the introduction of this space and its properties (see also Benth and Kriihner (9) for 
a further analysis of this space). 

Consider X defined in (13. Il l for H = H w and A = d/dx, the derivative operator. Then X can be 
considered as the dynamics of the forward curve, that is, f(t, x) := 6 x (X(t)) = X(t)(x), where f(t, x) := 
F(t, t + x), and t i-» F(t, T), t < T is the arbitrage-free forward price dynamics of a contract delivering 
an asset (commodity or stock) at time T (see Benth and Kriihner J3). We note that the semigroup of A 
will be the right shift operator S(t)f = /(• + t), and that 

S x S(t)g = g(t + x) = S x+t g. 

for any g £ H w . We find from the mild solution of X in (13.3b that 

(4.3) f(t,x) = f 0 (t + x) +S X I S(t - s)y 1/2 (s)dB(s) 

Jo 

where /q(£ + x) = S x S(t)Xo. Note that by Lemma 3.2 in Benth and Kriihner J3, it holds 

lim (X 0 ,S*(t)h x ) H = lim S x S(t)X 0 = lim f 0 (t + x) = / 0 (oo). 

t—> oo t—> oo t—t oo 

Here, fo(oo) denotes the limit of fo(y) as y -A oo, which exists. Hence, from the mild solution in (14.3b . 
the mean of /(£, x) for given x £ R+ has a limit /o(oo) as time tends to infinity. 

Now we investigate the stochastic integral in (14.3b in more detail. By Thm. 2.1. in Benth and Kriihner 13, 
there exists a real-valued Brownian motion b x such that 

(4.4) S x [ S(t - s)y 1/2 (s)dB(s) = f a(t,s,x)db x (s), 

Jo Jo 

where 

er 2 (f, s, x) = (■ d x S(t - s)y 1 ^ 2 {s)Q(S(t - s)3t 1 / 2 ( s ))*5*)(l) 

= S x S(t - s)(y 1 / 2 (s)Qy 1 / 2 (s))(6 x S(t - s))*(l) 

= 5 x+t _ s (y 1/2 Qy 1/2 (s))s* +t _ s (i). 

We remark that the Brownian motion b x depends on x, since the representation in Thm. 2.1. in Benth and 
Kriihner |9) is for a given linear functional, which in this case S x . We know that S*(l) = h x (-) (see e.g. 
Benth and Kriihner S3), and therefore 

(4.5) cr 2 (t, s, x) = {y 1/2 {s)Qy 1/2 (s)(h x+t s{-)))(x + t-s). 

Hence, we map the function h x +t-s by the operator y 1 ^ 2 {s)Qy 1 ^ 2 {s), and evaluate the resulting function 
in H w at x + t — s. As }4/ 2 is stochastic, we get a stochastic volatility cr(t , s, x), which is depending 
on current time t, previous times s and the ’’spatial” variable x. In particular, the spot price dynamics 
S(t) := f(t, 0) becomes 

S(t) = fo(t)+ I cr(t, s, 0) dbo(s). 

Jo 

I.e., the spot price dynamics follows a Volterra process where the integrand a(t. s, 0) is stochastic. We refer 
to Barndorff-Nielsen, Benth and Veraart 10 for an application to Volterra processes (and more specifically. 
Brownian and Levy semistationary processes) to model spot prices in energy markets. 

Let us carry our discussion further, and suppose that Q commutes with 4 'q and C(t) for t> 0, as well as 
that we have C(T)Q = C(T Q) and QC(T) = C(QT) for any T £ Ji. Then we recall from the previous 
Section that y(s) will commute with Q for every s > 0. The process 3 / ’ l,/2 (-s) will also commute with Q, 
and 

o 2 (t,s,x) = (y(s)Qh x +t- s )(x + t — s). 

Recalling the definition of 4 q in ( 13.7b . we find 

cr 2 {t, s,x) = 5 x +t-s(yQ(s){h x + t - s )) 

— Q,{,sJ)h'x-\-t—S’)h j x-\-t—s)w 

= (3^q(5), h x -\-t—s ® hx-\-t—s)l-L 
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since S z (f) = ( f,h z ) w for any / £ H w . Similar as in Prop. 12.51 we can calculate the cumulant of 
the process 4 ’q for any T £ H, and in particular we can calculate the cumulant of the process s (->• 
(yQ(s),h x+t -s <8> h x+t -s)-H for s < t by choosing T = h x+t - s ® h x+t - s ■ A simple calculation using 
the definition of h x in (14.2b shows that 

( h x +t—s ® hx+t—s )(/)=(/(° ) + J o f{y)dy S jh x-\-t—S — -£x-\-t—s (f)h x , 

where X x £ ff* is defined as l x {f) = 6 0 (f) + f'(y) dy for any f £ H w . 

In the above considerations we obtain a ’’marginal” dynamics, in the sense of a dynamics for a forward 
contract with fixed time to maturity x. We now represent the forward price dynamics as a space-time 
random field to emphasize also its spatial dynamics (i.e., its dynamics in time-to-maturity x). First, from 
Prop. 3.6 in Benth and Kriihner lUOl we find for any / £ H w , 

{s{t - s)y 1/2 (s)f)(x) = ( s(t - s)y 1/2 (s)f , h x ) w 

= (/, (S(t^s)y 1 / 2 (s))*(h x )) w 

= (s(t- s )y 1 / 2 (s)Y(h x mf(o) 

POO 

+ / w{y){S{t-s)y 1/2 {s)y{h x )'{y)f\y)dy 

J 0 

= (y 1/2 (s)s*(t-s)h x )(o)f(o) 

poo 

+ / w{y){y 1,2 {s)S*{t-s)h x )'{y)f\y)dy . 

J 0 

Again from Prop. 3.6 in Benth and Kriihner flO) , 

pOO 

S*(t)h x (-) = h x (0)(S(t)h.)(0) + / w{y)(S(t)h.)'(y)ti x (y) dy . 

Jo 

But S(t)h x (y) = h x (y + t) and h' x {y) = w~ 1 {y)l(y < x ). Hence, 

S*(t)h x (-) = h.(t) + [ w~ 1 {y + t)l{y + t < ■) dy 

Jo 

= /l ‘(') + It <-) d y 

= ^-£+;c(‘) • 

If we use the notation that B(ds , dy) := d x B(ds , y ) dy , we find 

s x [ S(t-s)y 1 / 2 (s)dB(s)= [ (^ 1/2 (s)/i x +t- s )(0)dH( S ,0) 

Jo Jo 

n oo 

w(y)(y 1/2 {s)h x+t - s )'{y ) B(ds , dy), 

or, a representation of /(f, x) := 5 x (X(t)) as a spatio-temporal random field 

pt pt pOO 

f{t,x) = fo(t + x)+ (y 1/2 {s)h x+ t- s )(0)dB(s,0) + / w(y)(y 1/2 (s)h x+t -sY(y) B(ds, dy). 

Jo Jo Jo 

Note that B(t , 0) = SoB(t) = (B(t), ho) w = ( B(t ), 1)^ is a real-valued Brownian motion with variance 
\Q 1/2 l\ 2 w . Hence, b 0 (t) := B(t, 0)/|Q 1/2 l|„, is a real-valued standard Brownian motion and we can view 
the first integral as an Ito integral of a volatility process given by s >-» (3^ 1 / 2 (s)/i a; +t_ s )(0) for s < t where 
a; is a parameter. It becomes a real-valued Volterra process with parameter x. For the second integral, we 
integrate with respect to a spatio-temporal random field (s,y) e -> B(s,y) over [0, t] x R + , thus becoming a 
stochastic Volterra random field. This part is analogous to an ambit field, a class of spatio-temporal random 
fields defined in Barndorff-Nielsen and Schmiegel 0. In a special case, the ambit fields take the form 

n oo 

g(t,s,x,y)rj(s,y)B(ds,dy) 
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for a stochastic random field 77 and a deterministic kernel function g. Under appropriate integrability 
conditions, the ambit field A(t, x ) is well-defined (see e.g. Barndorff-Nielsen, Benth and Veraart JT|). We 
observe that we can identify w(y) with the kernel function g, giving a very simple kernel. On the other 
hand, the volatility field 77 is more complex in /, as it is also ./'-dependent and not only s and y dependent. 
Our stochastic volatility model serves as a motivation for an extension of the ambit field models. We refer 
to Barndorff-Nielsen, Benth and Veraart El and 0 for an application of ambit fields to energy forward 
price modelling. 

We finally remark that in many commodity markets one observes an increasing volatility with decreasing 
time to delivery, known as the Samuelson effect (see Samuelson 11231 ). To include this in our dynamics of 
X, we can add an operator T(f) £ TL, possibly time-dependent, such that 

d,X(t) = AX(t) dt + dB(t ). 

Much of the analysis above can, under natural integrability conditions on T, be carried through for this 
model. 


Appendix A. A result on symmetric Hilbert-Schmidt operators 

In this section we provide the arguments for a claim in the proof of Proposition l2.10l for the convenience 
of the reader. 


Lemma A.l. Let H be a separable Hilbert space and let (ek)k£N an orthonormal basis of H. Then any 
symmetric Hilbert-Schmidt operator V on H can be written as 

V = ^2 7k,iek ® ei, 

k,ie n 

with Jf kl 7 ^ < 00 and 7 k ,i = 7 i,k- 

Proof. Recall that the space of Hilbert-Schmidt operators on H, denoted by %, is isometrically isomorph 
to H* ® //, which we identify as H ® H. Therefore, any Hilbert-Schmidt operator V can be written as 


v= £ 7/c,ze fc (g) ej, 

k,ie N 

for a sequence of constants ( r yk,l)k,le^- 

As for the square-summability of the constants, we have necessarily 


= y, ^2 Tk,i{ek ® ei)ei 

ign k,ien H teN isn 


H 


iGN ' i£N mGN 

= 'y ^ e m)H 

mGN 

= V 7 2 

/ j h,m 

i. mGN 


H 


So, in order for the norm to be finite, the double sequence {"/k,i)k,ieN has to be square-summable. 

As for the property of being symmetric, we need to have for all /, g £ dom(V) that 

(A.l) (Vf,g) H = (f,Vg) H . 

Since Hilbert-Schmidt operators are bounded (even compact), one has dom(V) = H. As a side-remark, 
this furthermore implies that dom(V*) D dom(V) = H, which in turn implies that symmetric Hilbert- 
Schmidt operators are already self-adjoint. The terms in (IA.lb can be evaluated to be 

(Vf,g) H =Y J 1k,ifkgi and (f, Vg) H = ^ lk,ifi9k , 
k,i k,i 
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where /*. = (/, , and similarly for gi. These terms can only be equal for all f,g£ dom(V) = H if 

either k = 1 or if 7 k,i = 7 i t k, which implies the assertion. □ 
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